perm filename QFACE.SAI[GEO,BGB] blob sn#083868 filedate 1974-01-26 generic text, type C, neo UTF8
COMMENT ⊗   VALID 00002 PAGES
C REC  PAGE   DESCRIPTION
C00001 00001
C00002 00002	BEGIN "QFACE"
C00005 ENDMK
C⊗;
BEGIN "QFACE"
	REAL AA,BB,CC,KK,QQ;
	REAL X0,Y0,Z0;
	REAL X1,Y1,Z1;
	REAL X2,Y2,Z2;
	REAL X3,Y3,Z3;

	X1←+1;	Y1← 0;	Z1←0;
	X2←-1;	Y2←+1;	Z2←0;
	X3←-1;	Y3←-1;	Z3←0;

	KK ← X1*(Z2*Y3-Y2*Z3) + Y1*(X2*Z3-Z2*X3) + Z1*(Y2*X3-X2*Y3);
	AA ← (Z1*(Y2-Y3) + Z2*(Y3-Y1) + Z3*(Y1-Y2));
	BB ← (X1*(Z2-Z3) + X2*(Z3-Z1) + X3*(Z1-Z2));
	CC ← (X1*(Y3-Y2) + X2*(Y1-Y3) + X3*(Y2-Y1));
	QQ ← 1/SQRT(AA*AA+BB*BB+CC*CC);
	AA ← AA*QQ;BB ← BB*QQ;CC ← CC*QQ;KK ← KK*QQ;

	OUTSTR(CVG(AA)&9&CVG(BB)&9&CVG(CC)&9&CVG(KK)&13&10);
	OUTSTR(CVG(CC*1 - KK)&13&10);
END "QFACE";